A Phase Front Instability in Periodically Forced Oscillatory Systems 
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Multiplicity of phase states within frequency locked bands in periodically forced oscillatory systems 
may give rise to front structures separating states with different phases. A new front instability is 
found within bands where u forcing /(^system — 2n (n > 1). Stationary fronts shifting the oscillation 
phase by n lose stability below a critical forcing strength and decompose into n traveling fronts 
each shifting the phase by n/n. The instability designates a transition from stationary two-phase 
patterns to traveling n-phase patterns. 



Periodic forcing of a single oscillator can lead to rich 
dynamics including quasiperiodic oscillations, frequency- 
locked bands ordered through the Farey construction, 
and low-dimensional chaos ■ A typical feature of a 
periodically forced oscillatory system is the multiplicity 
of phase states within a given frequency-locked band || . 
This feature becomes particularly significant in spatially 
extended systems where phase fronts separating different 
phase states may appear. The simplest situation arises 
in a system that is forced at twice the natural oscilla- 
tion frequency (hereafter the 2:1 band). A phase front 
(kink) connecting two uniform states whose phases of os- 
cillations differ by it then exists (hereafter a "7r-front"). 
The stability properties of this type of front are well 
studied As the strength of forcing is decreased 

a stationary front loses stability to a pair of counter- 
propagating fronts through a pitchfork bifurcation. The 
instability, known also as the Nonequilibrium Ising-Bloch 
bifurcation , is responsible for the destabilization of 
standing waves and the onset of traveling wave phenom- 
ena including spiral waves. 

The low resonance bands, 2:1 and 3:1, have been stud- 
ied both theoretically 111^] and experimentally p3| . 



All phase-front solutions in these bands shift the oscilla- 
tion phase by the same angle (in absolute value): n in the 
2:1 band and 2ir/3 in the 3:1 band. At higher resonance 
bands phase fronts that shift the phase by different an- 
gles may coexist, for example, 7r-fronts and 7r/2-fronts in 
the 4:1 band. In this Letter we report on a new front in- 
stability: Upon decreasing the forcing strength a 7r-front 
within the 2n:l band (n > 1) loses stability and decom- 
poses into n interacting n/n- fronts. We analyze in detail 
the 4:1 resonance case and bring numerical evidence for 
the existence of this type of instability in higher reso- 
nances. 

We consider an extended system that is close to a Hopf 
bifurcation and externally forced with a frequency about 



four times larger than the Hopf frequency. The amplitude 
of oscillations satisfies the parametrically forced complex 
Ginzburg-Landau (CGL) equation [ ]l4] , |l5| 

A t = (p + iv)A + (1 + ia)A xx - (1 - i[3)\A\ 2 A (1) 

where the subscripts t and x denote partial derivatives 
with respect to time and space, A(x, t) is a complex field 
and v, a, (3 and 74 are real parameters. We first consider 
the gradient version of Eqn. ([!]) by setting v = a = (3 — 
and then rescale time, space and the amplitude A by [a, 
\J fJb/2 and l/^/Jl, respectively. Keeping the same nota- 
tions for the scaled quantities the gradient version takes 
the form 



A t = A+ -A xx - \A\ 2 A + 74A* 



(2) 



For < 74 < 1 Eqn. (g) has four stable phase states: 
A±i = ±A and A±i — ±«A, where A = — 74. Front 

solutions connecting pairs of these states divide into two 
groups: 7r- fronts connecting states with a phase shift of 



= A + \ tanhx , 
A-i^+i = A + i tanhcc , 



(3) 
(4) 



and 7r/2-fronts connecting states with a phase shift of 7r/2 
(see Fig. ||a). When 74 = 1/3 the 7r/2-fronts are given 

by 

A+i^+i = \ + i 1 ~ «)tanha;] , 

A-i^+i = i + (l + i)tanha;] , (5) 

A-i^-i = -A + i^ +i and A +i ^._i = -A-i-> +1 . Ad- 
ditional front solutions follow from the invariance of 
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Eqn. (|2|) under reflection, x — > — x. For example, the 
symmetric counterparts of A + i^ + \(x) and A + \^-i(x) 
are A +1 ^ +l (x) = A +i ^ +1 (—x) and (x) = 

A+i^-^-a;). 

Stability analysis of the 7r-fronts (|J) and (|4|) shows that 
they are stable for 74 > 1/3. To study the instability at 
74 = 1/3 we rewrite Eqn. (|) in terms of U = $t(A)+%(A) 
and V = SR(A) - 3(A) : 

tf t = ^ + ^-~tf 3 -~(tf 2 -3F 2 )E/, 

^ = ^+^-^ 3 -^ 2 -3t/ 2 )^, (6) 
where 

d = 74 - 1/3 . 

At d — the two equations decouple and assume the so- 
lutions U = <j\Aq(x — x\) and V — a 2 Ao(x — x 2 ), where 

Aq = y^t&rihx, o\ 2 = il an d X\ and x 2 are arbitrary 
constants. Consider now d 7^ but small. The coupling 
between U and makes x\ and 2:2 slow dynamical vari- 
ables and U and V can be written as 

U = <tiA [x — Xi{t)} + u , 

V = a 2 A [x- x 2 (t)] +v, (7) 

where u and v are corrections of order d. Inserting these 
forms in Eqns. (|^) we obtain 

Ji\u = <JiXiA' (x — xi) 

- -dcri - xi) - 3Al(x - x 2 )]A (x - xi) , 

where 7ii = — 1 — |^t + 2Aq(x— Xi). A similar equation 
is obtained for u with the indices 1 and 2 interchanged. 
Solvability conditions lead to the equation 

Xi = d [ dx 

16 J-oo 

x tanh(x — xi) sech 2 (x — Xi) tanh 2 (x — x 2 ) , (8) 

and to a similar equation for x 2 with the indices 1 and 2 
interchanged. Defining a translational degree of freedom 
C = |(xi + X2) and an order parameter \ — 5(^2 ~ x i)i 
we obtain from (@) 

C = 0, x = -g<iJ( x ), (9) 

where 

J( x ) = I( a ) = 6(a~ 1 - a' 3 ) + (1 - 3a~ 2 )G(a) , 
G(o) = (1- a- 2 ) In(^), 

and a = tanh2%- Note that Eqns. (^) are valid to all 
orders in x an d to linear order around 74 = 1/3. 



Figure || shows the potential V(x) — j X I{z)dz 

associated with Eqn. (^) for d > (74 > 1/3) and d < 0. 
There is only one extremum point, \ — 0, of V. For 
d > it is a minimum and x converges to zero. For 
d < it is a maximum and x diverges to ±00. To see 
how the x dynamics affect the front solutions we rewrite 
the solution form (^) in terms of the amplitude A and 
bring it to the form 

A(x,t) = A^i^ + i(x - xi) + A +1 ^ +i (x - x 2 ) - X + R, 

(10) 

where A-i->+x and A + \^ +i are given in Eqns. (||), R 
is a correction term of order d (related to u and v in 
Eqns. (]?])) and for concreteness we chose o\ = —o 2 = 1. 
Equation ([n]) describes two interacting 7r/2-fronts cen- 
tered at x\ and x 2 . When d > \x\ decreases in 
time and the two 7r/2-fronts attract. As \ — ► (or 
x 2 — ► xi) the two 7r/2-fronts collapse into a single n- 
front, A(x, t) = A-i^ +i (x — x\) + R, given by Eqn. (^). 
When d < \x\ increases in time and the two 7r/2-fronts 
repel. Perturbing the unstable x = solution, the sin- 
gle 7r- front decomposes into a pair of 7r/2- fronts. Fig. ||a 
shows phase portraits of the tt and 7r/2-fronts (dashed 
lines) and the time evolution of an unstable 7r-front for 
d < (solid lines) obtained by numerical integration of 
Eqn. (||) . The approach of the phase portrait to the fixed 
point A + i on the $s(A) axis describes the decomposition 
into a pair of 7r/2-fronts. This behavior persists arbi- 
trarily close to d = and is related to the absence of 
minima in the potential V(x) for d < (see Fig. |l|fe). 
At d — there exists a continuous family of stationary 
solutions describing frozen (non- interacting) pairs oiw/2- 
fronts with arbitrary separations X2 — x\. This solution 
family spans the whole phase space inside the dashed tri- 
angle in Figure |^a. Because of the parity breaking sym- 
metry x ~ * ~X each 7r-front may decompose into one of 
two pairs of it/2- fronts with phase portraits approaching 
the fixed points A + i and 

We extended the derivation of Eqns. (j|) to the non- 
gradient equation (Q), treating the constants v, a, /3 as 
small parameters. The x equation remains unchanged. 
The C equation takes the form 

— C = vF„(x)+aF a ( X )+PF p (x), (11) 
0-1 

where 

F a = \l{a), 

F, = 3a-^l- Z -a- 2 )- 9 -a-"G(a). 

Notice that F u , F a , and Fp are odd functions of x an d d° 
not vanish when d = 0. When \x\ — > 00 the right hand 
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FIG. 1. The potential V(x)- (a) For d > the extremum 
at x = is a minimum and and x converges to 0. (b) For 
d < the extremum is a maximum and x diverges to ±oo. 



side of (fLlJ) converges to i|(z/+/3), the speed of a 7r/2-front 
solution of Eqn. (Q). The x = solution (representing 
a 7r-front) remains stationary (£ = 0) in the nongradient 
case as well. At 74 = 1/3 (d = 0) it loses stability and 
decomposes into a pair of 7r/2-fronts which approach the 
asymptotic speed |(z/ + /3). Depending on the initial sign 
of x the pair may propagate to the left or to the right. 

This behavior is different from that near the Noncqui- 
librium Ising-Bloch front bifurcation within the 2:1 band. 
In that case, a stationary Ising front loses stability to a 
pair of counter-propagating Bloch fronts in a pitchfork 
bifurcation. Associated with the bifurcation is a transi- 
tion from a single- well potential (Ising front) to a double- 
well potential (pair of Bloch fronts). A comparison with 
the potentials in Fig. [I] shows the essential difference be- 
tween the two front instabilities. In the 2:1 band the 
Bloch fronts approach the Ising front and coincide with 
it as the distance to the bifurcation point diminishes to 
zero. In the 4:1 band, on the other hand, the asymptotic 
solutions just below 74 = 1/3 (the 7r/2-front pairs as 
|x| — ► 00) are not smooth continuations of the stationary 
7r-front (the x — solution) at 74 = 1/3. In particular 
their speed remains finite (%{y-\-0)) as 74 approaches 1/3 
from below. At 74 = 1/3, a whole family of propagating 
solutions appears with speeds ranging continuously from 
+ /?) to zero (pertaining to 7r/2-front pairs separated 
by distances ranging from infinity to zero). 

The instability of n- fronts at 74 = 1/3 (d = 0) de- 
termines the structure of stable periodic patterns below 
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FIG. 2. Decomposition of n- fronts into n/n- fronts within 
2n:l bands for (a) n — 2, (b) n — 3, and (c) n = 4. The 
dots are the uniform phase states along the circle of constant 
amplitude \A\. The 7r-fronts are the dashed lines connecting 
the states A+i and A-i on the K(A) axes. The 7r/ra-fronts 
are the dashed lines connecting the adjacent points along the 
circle. The time evolution of the front decomposition obtained 
by numerical integration of Eqn. ( p"2| ) is shown as the series 
of solid lines. Parameters: all undetermined coefficients in 
Eqn. ( |l2] ) were set to zero except as indicated below, (a) 
74 = 0.3, ^4 = -1.0, (b) 7 6 = 0.9, m = -1.0, Me = -1.0, (c) 
78 = 0.75, fi4 — —0.5, fj,6 — —0.5, fis = —1.0. 
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FIG. 3. Collapse of a rotating four-phase spiral- wave into a 
stationary two-phase pattern. The left column is and the 
right column arg(A). (a) The initial four-phase spiral wave 
(computed with 74 < 1/3). (b) The spiral core, a 4-point ver- 
tex, splits into two 3-point vertices connected by a tt- front, 
(c) A two-phase pattern develops as the 3-point vertices fur- 
ther separate, (d) The final stationary two-phase pattern. 
Parameters: 74 = 0.6, v = 0.1, x = [0, 64], y = [0, 64]. 



and above the instability. In the range 74 > 1/3 two- 
phase patterns, involving domains separated by 7r-fronts, 
prevail. Below 74 = 1/3 four-phase patterns dominate. 
Four-phase patterns are not stable for 74 > 1/3 despite 
the stability of the 7r/2-fronts because of the attractive 
interactions among these fronts. In the gradient case 
[Eqn. (||)] all solutions are stationary. In the nongradi- 
ent case [Eqn. (|l|)] the two-phase patterns are stationary 
while the four-phase patterns propagate. Fig. ^p, shows 
a grey-scale map of a rotating four-phase spiral wave for 
74 < 1/3. Figs. |]&, c, d show the collapse of this spi- 
ral wave into a stationary two-phase pattern as 74 is in- 
creased past 1/3. The collapse begins at the spiral core 
where the 7r/2-front interactions are the strongest. As 
pairs of 7r/2-fronts attract and collapse into tt- fronts, the 
core splits into two vertices that propagate away from 
each other leaving behind a two-phase pattern. 

To test whether the instability of 7r-fronts exists at 
higher resonances we integrated numerically the higher 
order equation (!(]] 

A t = l -A xx + A + ^\A\ 2 A + fi^A + v 8 \A\ 6 A (12) 
+ 7 4 A* 3 + 7 6 A* 5 + 78 A* 7 . 

Indeed 7r-fronts within the 6:1 and 8:1 bands become 
unstable as 7 6 and 73, respectively, are decreased be- 
low a critical value. Fig. ||t> shows the decomposition 
of a 7r-front into three 7r/3-fronts within the 6:1 band, 
and Fig. shows the decomposition into four 7r/4-fronts 
within the 8:1 band. Our conjecture is that the instabil- 
ity is general, occurring within any 2n:l band (n > 1). 

The phase front instability and the associated tran- 
sition from stationary two-phase patterns to traveling 
four-phase patterns within the 4:1 band may be tested 
in experiments on the ruthenium-catalyzed Belousov- 
Zhabotinsky reaction subjected to periodic (in time) uni- 
form illuminations fl3|| . 
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